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Abstract: The superconducting transition is studied within the one-loop RG in fixed di- 
mension D = 3 and at the critical point. A tricritical behavior is found and, for k > k c , 
an attractive charged fixed point, distinct from that of a neutral superfluid. The critical 
exponents of the continuous transition are evaluated and it is shown that the anomalous 
dimension of the gauge-field equals unity. This implies the proportionality of the magnetic 
field penetration depth and the superconducting correlation length below the transition. The 
penetration depth exponent is non-classical. We argue that it cannot be extracted from the 
dual theory in a straigthforward manner since it is not renormalized by fluctuations of the 
dual field. 
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The problem of a charged scalar field coupled to a gauge vector potential arises frequently 
in theoretical physics. In its original version, it describes formation of Meissner state in 
superconductors |] and Higgs mechanism in particle physics [[J. Furthermore, the nematic- 
smectic-A transition in liquid crystals |§ and, more recently, the transitions between plateaus 
in the quantum Hall effect @, H and the finite-magnetic-field critical behavior in extreme 
type-II superconductors have also been related to this problem. In a superconductor, 
the scalar field represents the fluctuating superconducting order parameter \& which, BCS 
pairs being charged, is coupled to fluctuations in the electromagnetic potential A. At the 
mean-field level the transition is discontinuous and remains so when fluctuations in \1/ are 
included via the e(= 4 — Z))-expansion [jlj. Numerical simulations 0, || of related lattice 
models support this scenario for smaller values of the Ginzburg parameter k. However, for 
large k, the results are consistent with a continuous, second-order phase transition [^, [|. The 
picture obtained in numerical work is in accordance with the l/n-expansion |1|, [ij and with 
duality arguments which connect the lattice version of the theory to a dual gas of interacting 
vortex loops and the "inverted" 3D XY model |7], [K| . 

In this Letter we study this superconducting-Higgs electrodynamics (SHE) directly in 
D = 3 within perturbation theory at the critical point corresponding to the charged super- 
fluid. Our results are as follows: We first show that the anomalous dimension of the gauge 
field, t]a, equals unity to all orders in perturbative expansion. By combining this result 
with the Josephson relation, we argue that the magnetic field penetration depth, A, and the 
superconducting correlation length, £, diverge with the same exponent (v) as the transition 
is approached from below. These results should be contrasted with t]a = and A oc a/£ a t 
the unstable critical point for neutral superfluid. We then demonstrate that our one-loop 
results imply the presence of the tricritical point in SHE. For bare k < k c the renormalization 
group (RG) flows are unstable and the transition is likely first order while, for k > k c , we 
find a stable fixed point indicating a continuous transition. By selecting the regularization 
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which reproduces the established numerical value of k c , we evaluate v and the anomalous 
dimension of 77^ . Our results strongly suggest that the exponent for A is non-trivial and 
therefore different from the mean-field value suggested by Kleinert et al. |10| We propose 



that this difference has a physical origin and that this form of dual theory does not offer 
any simple way of determining the exponent of A. Since the penetration depth is of direct 
experimental interest [|11], |T2|| , the description of the critical behavior within the original 
SHE theory retains its physical significance. 

We are interested in the infrared behavior of the SHE: 

H = |^|(V-zeA)^(f)| 2 + /i 2 |^(r)| 2 + ^(r)| 4 + i(Vxl) 2 ] , (1) 

where /i 2 oc (T — T c0 ), T c0 is the mean-field transition temperature, b is a temperature 
independent constant, and e is the charge of a BCS pair. For generality, we assume that the 
order parameter \1/ has n complex components and that the system is D-dimensional, n = 1 
and D = 3 being eventually the case of physical interest. We chose to work in the gauge 
where the vector potential is purely transverse, i. e. where the bare gauge-field propagator 
is: Dij(q) = (5ij — qi(jj)/q 2 - First, let us discuss the non-perturbative results concerning the 
anomalous dimension of the gauge-field and the divergence of the penetration depth. The 
anomalous dimension of the gauge-field propagator is defined as: 

d\ogZ A 

Va = ~ hm — : — pr- 2 
p^o dlogyo) 

where Za is the gauge-field renormalization factor and p is the momentum of the gauge-field 
propagator. The /5-function for the charge is 

(3 = d % =c 2 (D 1 dl °S ZA \ (3) 

where e 2 = e 2 /p is the dimensionless renormalized charge. On approaching the critical point, 
e r — > eo and (3 e —>■ 0, so taking the limit p — > in Eq. (3) we obtain: 

e 2 (D- 4 + VA ) = . (4) 



Assuming a stable charged fixed point in the theory eo 7^ one obtains 

VA = A-D . (5) 

This exact result has an important physical consequence, since it modifies the relation be- 
tween the correlation length and the penetration depth as the critical point is approached 
from below. In general, the penetration depth scales with the superfluid density below the 
transition as: 

A^- 2 ocp s , (6) 

where p s oc ^ 2 ~ D is the Josephson relation Close to the transition controlled by the 

attractive charged fixed point it follows that 

A oc £ (7) 

for all D. In contrast to the scaling governed by the XY fixed point, where t]a = and 
A oc £( D - 2 )/ 2 ; the ratio between the two lengths close to the charged critical point approaches 
a finite constant. The divergences of both lengths are determined by the same exponent v. 

To obtain the flow-diagram for the coupling constants and the values of critical ex- 
ponents one must rely on some approximation for the /3-functions. Here we perform the 
perturbative calculation of /5-functions in fixed dimension and at the critical point. To the 
lowest order in perturbation theory, the contributions to the self-energy, polarization and 
the quartic vertex are given by the diagrams in Fig. 1. Due to the choice of gauge, the 
remaining one-loop diagrams for the quartic vertex all vanish when the external momenta 
go to zero 0. Note that this procedure explicitly preserves Ward identities associated with 
gauge invariance. Since we wish to work directly in D = 3, we are forced to define the 
renormalized value of b at a finite momenta of external legs to avoid the infrared divergence 
in the last diagram on Fig. 1. This divergence is a consequence of gauge invariance which 
requires massless gauge-field. The renormalized coupling constant b r is defined at the usual 
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symmetric point: 



h-ki = ^\^p\ i,j = 1,2,3 (8) 



and at the critical point where the renormalized mass of ^ vanishes. There are two relevant 
coupling constants in the problem: the quartic term coupling and the charge. The standard 



procedure [14] gives the renormalized coupling constants to lowest order: 



, , , 2nT(l-D/2)T 2 (D/2) . „_ 4 



e 2 = e 2 + v 1 ' K ' 'e k a u ~* (9) 



n + 3 F(2-D/2)r 2 (D/2-l) D _ 4 



v 2 2 -^/ 2 7 (4n) D / 2 T(D - 2) 
4QD - l)r(2 - £>/2)r(D/2 - ljTQD/g) „_ 4 _ 2(D - l)r(2 - D/2)T 2 (D/2 - 1) D _ 4 
(4tt) d / 2 T(D - 1) P (47r) D / 2 r(D - 2)2 2 - D / 2 P ' 

where T(x) is the factorial function, the momentum scale p is defined above and q is the 
momentum of the gauge-field propagator. Hereafter, we set D = 3 and n — 1 and define 
dimensionless couplings with respect to the external momentum p: b = b/p, e 2 = e 2 /p. If 
we choose the reference momentum of the gauge-field propagator as q = p/c, where c is a 
constant, the /3-functions are: 

A ^^=- s+ ^ +1 ^ 2 -^ 2+ 57f 4 • (11 > 

and we dropped the subscript "r" for renormalized coupling constants in the last two equa- 
tions. 

The above /^-functions explicitly depend on the introduced ratio of momenta, c. This 
reflects the known property of the RG in fixed dimension that, unlike in the e-expansion 
scheme, the coefficients in the perturbative series for /3-functions are procedure dependent 
T4j . Consequently, an additional information is needed to fix the value of c in the one-loop 



calculation. We now note an important property of Eqs. (10,11): If c is treated as a free 



parameter, the RG flow diagram has a structure pictured in Fig. 2 (for c > 5.16). Besides 
the standard Gaussian and the neutral superfluid fixed points, there are two charged fixed 
points of the above flow equations: (&_,eg) and (fo + ,e§), where <3q = 16/c and i> + and 6_ 
are the real roots of the equation /?&(&, eo) = 0. Stability analysis shows that the fixed point 
with the larger value of b = b + is attractive, while the one with 6 = 6_ is unstable in the 
direction of quartic term coupling. Now we look for the straight RG trajectories, b = 2/t 2 e 2 , 
by demanding the invariance of the constant k under RG transformation: 

d(b/e 2 ) 



d\og(p) 

In our one-loop analysis this leads to the equation: 



. (12) 



ft(2«^,e*) = , (13) 



resulting in two straight-line RG trajectories for which 



<- = W ' (14) 



Together with the result for the stability of the (6_,Cq) fixed point, this implies that there 
is a tricritical line in the theory given by a GL parameter: 



, L c + 8 -yjc 2 + 16c -32(2 + ^2) 



kT = —^ = 1 t= . (15) 

c 2eg 8(2^2 + 1) 1 ^ 



This solution appears physically plausible and is in accord with Ref. ]TB| in that the tricritical 
point is determined by a particular value of the Ginzburg parameter, as one would expect 
from a mean-field argument. Consequently, we fix the value of c by demanding that k c in 



Eq. (15) matches the value obtained via the duality transformation of the lattice SHE [15| 
(see also numerical Monte Carlo results of Ref. f|): k c = 0.8/a/2. This condition results 
in c = 5.7 (> 5.16). The other straight trajectory given by the k+ solution connects the 
Gaussian and the superconducting fixed point. 
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Having fixed the value of the parameter c we thus obtain the flow diagram of Fig. 2 
16fl , [[HJ. The exponents at the attractive fixed point are 

77* = ^ = -0.70 , (16) 

- = ^(l + ^"f) =0-53 , (17) 

with the numerical values calculated for c = 5.7, (k c = 0.8/^/2). It is worth mentioning that 
our procedure gives a respectable value for the correlation length exponent at the neutral 
superfluid fixed point: vxy = 0.63. The value of the anomalous dimension is rather large and 
negative but it does satisfy 77$ > — 1 in D — 3. Note that small reductions in the assumed 
value for n c rapidly make 77$ less negative (as do the next order perturbative terms) while the 
value for v is more robust. For example, k c = 0A5/\^2 || gives r]^, = —0.20 and v = 0.62. 
The other exponents follow from standard scaling relations. It is conceivable, however, that 
the hyperscaling relation does not hold, due to the presence of long-range gauge forces in (1). 
In that case we could define a characteristic dimension, d c , from 2 — a = d c v. By combining 
our results with the prediction of the dual theory a = axy = —0.013 one obtains d c = 3.8. 
This is close to d c = 4 which would arise from a naive scaling of current-current interactions. 

To the lowest order in e = 4 — D the /3-functions derived from Eqs. (9) completely 
reproduce the results of the RG defined at zero external momenta and finite mass |], [16| . 
This is to be expected since the dimensional regularization and the minimal subtraction 
scheme lead to unique values of the coefficients in the /^-functions ||14|| . The parameter c 
then does not appear at all, the attractive fixed point exists only for n > 182.95 and the 



exponents agree with ours in this limit. |L8 



The exponent for the magnetic field penetration depth is v = 0.53(0.62). While this 
appears close to the mean field value of 1/2 suggested by Kleinert et al. [l(| on the basis of the 
dual theory, it is clear from our procedure that this exponent is non-trivial and determined 
by the structure of the charged fixed point. This is an important aspect of the physics of this 
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problem: In the dual approach the partition function of the original problem is related to the 
one for the interacting vortex loops 0. The dual description is a theory of a scalar complex 
field coupled to a massive vector potential [0IJ and is in the universality class of the 3D XY 
model ||10|| . This description could be useful in providing the information on thermodynamic 
quantities. For example, the specific heat exponent a can be calculated from the hyperscaling 
relation, which holds in the dual theory, and has the value for the 3D XY model. The dual 
approach, however, offers little help in calculating exponents that characterize correlation 
functions of the original superconducting order parameter. In particular, the magnetic field 
penetration depth appears in the dual description via the mass of the "dual" vector field. If 
one assumes that this mass vanishes fast enough at the transition at the bare level, then, at 
the critical point of the dual theory, this mass does not renormalize, i.e. it remains equal 
to its bare value. [19j Since the mass of the dual vector field is what determines the inverse 
of the penetration depth, the implication is that the divergence of the penetration depth at 
the superconducting transition cannot be extracted from the RG analysis of the dual model, 
even in principle. To calculate the penetration depth exponent it is necessary to apply the 
RG analysis to the original SHE theory. This fact supplied the primary motivation for the 
present work. 

Our value for the penetration depth exponent, v = 0.53, is not seen in the experimental 



results of Ref. fllTl, which are consistent with the neutral 3D XY behavior, but is in excellent 



agreement with those of Ref. ||12|| . While the critical region in which charge of BCS pairs is 
relevant is rather narrow (even for HTS) due to the small value of the fine structure constant, 
the fact that our exponent is quite close to the mean-field value may effectively "enlarge" 
this region. In general, as we approach the transition from below, we expect that the mean- 
field behavior of the penetration depth crosses over to the intermediate regime of neutral 
superfluid followed by the ultimate charged superfluid critical behavior. The exponent in 
this intermediate 3D XY regime is ~ 1/3, considerably different from both the mean-field 
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value and the ultimate value of ~ 0.53. Consequently, on purely empirical grounds, the 
region over which the penetration depth can be described by the crossover behavior should 
be quite narrow. Finally, we should mention that experiments on other systems j|, [5], § 
might be even more promising in studying the SHE critical behavior since there the effective 
"fine structure constant" can be of order unity. 

In summary, we studied the superconducting transition by calculating the one-loop (3- 
functions in fixed dimension D = 3. The gauge invariance leads to infrared divergences which 
are handled by defining the renormalized coupling constants at finite values of the external 
momenta and right at the critical point. We eliminate the remaining freedom in the definition 
of renormalized couplings by requiring that our analysis yields to the previously established 



Tq , || numerical value of the Ginzburg parameter which characterizes the tricritical point. 



We then evaluate critical exponents at the attractive charged fixed point of our theory. It is 
shown that the penetration depth and the correlation length diverge in the same way close 
to the charged fixed point, to all orders in perturbation theory. |20| 
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Figure Captions: 

Figure 1. Lowest order contributions to the self-energy (a), polarization (b) and the 
quartic vertex (c). The full and dashed lines are the order parameter and the gauge- field 
propagators, respectively. 

Figure 2. The schematic flow-diagram for the dimensionless charge e and the dimen- 
sionless quartic term coupling b (for c > 5.16). Note that, for c = 5.7, = 2.81, b + = 3.22, 
6_ = 1.80, k + = 1.07/V2, K- = 0.80/ \/2. 
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